SU(3) Kondo effect in spinless triple quantum dots 
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We discuss a device — a purely capacitively coupled interacting spinless triple quantum dot system — for the 
observation of the SU(3) Kondo effect. Unlike more familiar SU(2) and SU(4) Kondo effects in quantum dot 
devices which lead to unitary linear conductance at low temperatures, the SU(3) Kondo scenario can be easily 
identified by the conductance pinned to a characteristic value of 3/4 of the unitary limit. This is associated 
with the interesting fact that the SU(3) Kondo effect does not occur at the particle-hole symmetric point, where 
the system is found instead in the valence-fluctuating regime with the total dot occupancy flipping between 
1 and 2, but for gate voltages in the two Kondo plateaux where the dot occupancy is pinned to an integer 
value, either 1 or 2. From the thermodynamic analysis in the Kondo regime we find that the effective impurity 
orbital moment, defined through the impurity orbital susceptibility (Xim P ) multiplied by the temperature, is 
T"ximp = 1 at high temperatures and then it increases to the characteristic value of Tximp = 4/3 corresponding 
to the three-fold degenerate local-moment fixed point where the impurity entropy is 5imp = In 3. Then, at 
much lower temperatures, the system flows to the non-degenerate strong-coupling fixed point in which the 
SU(3) Kondo effect takes place. We also report results about the robustness of the SU(3) Kondo effect against 
various perturbations present in real experimental setups, namely, unequal reservoir-dot tunneling couplings, 
gating effects and non-vanishing interdot tunneling rates. Finally, we describe possible mechanisms to restore 
the SU(3) Kondo physics by properly tuning the on-site dot potentials. We briefly comment on the spinfull 
case which has very different behavior and shows Kondo plateaus in conductance for all integer values of the 
occupancy, including at the particle-hole symmetric point. 

PACS numbers: 72.10.Fk, 72.15.Qm 



I. INTRODUCTION 

In metals, magnetic impurities are responsible for the 
anomalous behavior of the resistivity at low temperatures. 1 
Magnetic interactions result from high-order correlated tun- 
neling events of electrons that are hopping in and out of a lo- 
calized impurity. In this fashion, the impurity spin is screened 
through the formation of the Kondo spin singlet state. In semi- 
conductor quantum dots, the Kondo effect can also take place 
as theoretically predicted 2 -^ and experimentally observed?^- 
However, whereas the resistivity of a magnetically doped 
metal increases when the temperature is lowered below the 
Kondo temperature, Tk, in quantum dots (QDs) the linear 
conductance increases and eventually reaches its maximum 
value G = 2e 2 /h at zero temperature! 7 - This is due to the 
enhancement of the scattering rate which results in the open- 
ing of a channel for electrons that are perfectly transmit- 
ted through the QD. The main advantage of exploring the 
spin-1/2 Kondo effect in manufactured nanostructures such 
as QD a 5 ' 6 - 8,9 is their high tunability. Using gate electrodes 
it is possible to vary in a controllable manner the number of 
trapped electrons and the strength of tunnel coupling between 
the reservoirs and the localized dot states. Furthermore, as an 
additional advantage, we mention that QDs constitute perfect 
laboratories to test many-body effects under non-equilibrium 
conditions^— The influence of external fields such as finite 
bias voltages-l^ill or time varying ac fields 1 ' allows the obser- 
vation of the Kondo effect out of equilibrium. Remarkably, 
magnetic fields dramatically affect the Kondo state, even at 



equilibrium conditions. The level degeneracy, required for the 
formation of the usual spin-1/2 Kondo effect, is lifted in the 
presence of a magnetic field. There are, however, some ex- 
ceptions where a magnetic field facilitates the development of 
a Kondo state, for example in vertical QDs with an even num- 
ber of electrons. Here, the integer-spin Kondo effect occurs 
when the singlet and triplet states become degenerate because 
of the presence of the magnetic field 4^2°. Therefore, in gen- 
eral, magnetic fields either remove the Kondo effect or facili- 
tate its observation through the level degeneracy requirement. 

The fabrication of more complex nanostructures has opened 
new possibilities for the study of unconventional Kondo ef- 
fect. There exists a great variety of artificially fabricated sys- 
tems exhibiting exotic kinds of the Kondo effect: to men- 
tion just a few, nanodevices based on carbon materials like 
fullerenes, carbon nanotubes, and graphene^— Others sys- 
tems, such as p-doped nanowires, have demonstrated the 
Kondo effect assisted by holes^i The search for highly sym- 
metric Kondo singlets has been revived in carbon nanotubes 
and vertical double dot systems with the experimental demon- 
stration of the SU(4) Kondo effecti 23 - 25 ' 32 ~ ^l In carbon nan- 
otubes, the valley isospin together with the spin degree of free- 
dom manifests as a four-fold shell structure in the Coulomb 
blockade regime. 35 ' 36 In the low-temperature regime the fluc- 
tuations among the four quantum states lead to the observation 
of the SU(4) Kondo effect^ 3 - 3 ^ So far, the SU(2) and SU(4) 
Kondo effects have been extensively studied. There is, how- 
ever, very few works devoted to other possible symmetries for 
a Kondo singlet, and in particular to the SU(3) Kondo effect. 
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In a recent work 3 -^ the SU(3) Kondo effect has been suggested 
to be observable in triple quantum dot (TQD) in the quantum 
Hall regime. Here, we discuss a different setup — a spinless 
TQD with equal capacitative coupling Vi between all quantum 
dot pairs, as shown in the sketch in Fig.Q] as a suitable system 
for the observation of the SU(3) Kondo effect. The SU(3) 
Kondo physics takes place when there is a single electron or a 
single hole in the whole system, i.e., when either n = 1 (one 
electron and two holes) or n = 2 (two electrons and one hole). 
This defines three possible flavors corresponding to the posi- 
tion of the electron (or hole) in one of the three dots or leads. 
In this work, we will refer to the flavor degree of freedom also 
as the channel or orbital degree of freedom; for our spinless 
model, these expressions are fully interchangeable. Each dot 
is connected to two contacts in such a way that the tunnel- 
ing events conserve the flavor degree of freedom. Notice that 
in principle, this setup can be easily generalized to build an 
arbitrary SU(7V) Kondo state^d^ although that would entail 
designing a device with equal capacitive coupling between all 
QD pairs. Importantly, the only interaction among the dots is 
capacitative and there is no particle exchange from one dot to 
the others, i.e., the interdot tunneling is not allowed since this 
would destroy the flavor conservation rule. 

The recent progress in fabricating highly tunable TQDsii^ 
aims to provide a platform for testing a variety of predicted 
novel quantum information processing functionalities and 
many-body effects42i2i The goal of this work is to analyse the 
transport and thermodynamic properties of the SU(3) Kondo 
effect in highly symmetric capacitively-coupled TQD. Simi- 
lar study has very recently been performed in Ref. |5jJ; where 
comparison can be made, our results agree with theirs. We 
also study the effect of local perturbations on the SU(3) sin- 
glet Kondo state such as asymmetrical lead-dot couplings, 
finite interdot tunneling rates, non-equal charging energies, 
etc. Generally, these perturbations destroy the SU(3) singlet 
Kondo state, however, we propose a way to restore the Kondo 
resonance by properly gating the dot levels. 

In order to investigate in a general framework the different 
regimes encountered for the the TQD system we consider the 
operators of the SU(3) Lie algebra which describe the orbital 
(flavor) degree of freedom of the electrons. The thermody- 
namics analysis is performed by calculating the impurity or- 
bital (flavor) susceptibility Ximp{T), and the impurity entropy 
5'imp(7 1 )- Our results indicate that SU(3) Kondo physics oc- 
curs when the QDs are tuned to single occupancy, n = 1, 
or double occupancy, n = 2, which is achieved away from the 
particle-hole symmetric point, contrary to what happens in the 
more familiar SU(2) and SU(4) Kondo cases where at half 
filling the Kondo effect is present. Consequently, the Kondo 
peak itself is not symmetric with respect to the chemical po- 
tential as visible in the spectral densities for each dot. These 
values of the occupancy have important consequences for the 
linear conductance. In accordance with the Friedel-Langreth 
sum rule, the linear conductance is G = Go sin 2 6 where the 
scattering phase shift 6 is approximately given by 5 = irn/N; 
here n = (n) is the total TQD occupation and N = 3, while 
Go is defined as Go = e 2 /h (note that we are considering a 
spinless system, thus the spin factor 2 is not present in Go). 



Therefore, in the SU(3) Kondo regime with n = 1 and n = 2 
one has22 
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G K ondo = | e A" (!) 

This result must be compared to that at the p-h symmetric 
point where 5 = ir/2, and thu a 52 ' 53 

G p _ h = e 2 /h . (2) 

This paper is organized as follows. In Sec. I we introduce 
the model Hamiltonian to describe the TQD setup and discuss 
the theoretical tools to solve it. Section II is devoted to the 
study of the emergence of the SU(3) Kondo regime as a func- 
tion of various parameters, namely the dot level position (e^), 
the interdot Coulomb interaction (Vi), and the lead-dot tunnel- 
ing couplings (Vi). The discussion is based on the thermody- 
namics and we investigate the behavior of the impurity orbital 
susceptibility Xim P (T) and the impurity entropy S , ; mp (T). In 
Sec. Ill we list some signatures of the SU(3) Kondo state in 
the transport measurements. In Sec. IV we study the robust- 
ness of the SU(3) orbital Kondo singlet against diverse per- 
turbations, namely, asymmetric lead-dot tunneling couplings, 
different on-site energy values and possible leaking effects de- 
scribed by nonzero inter-dot tunnelling rates. In Sec. V we 
briefly consider the generalization to the spinfull problem and 
discuss the different kinds of the Kondo effect expected in that 
case. Finally our main conclusions are summarized in Sec. 
VI. 



II. MODEL AND METHODS 




Figure 1 : (Color online) Schematic representation of the capacitively 
coupled triple quantum dot system. Each quantum dot is attached to 
two electron reservoirs. We assume that a sufficiently large exter- 
nal magnetic field is applied to fully polarize the electrons, so that 
we may consider the system to be spinless. The only interaction 
between each pair of quantum dots is purely capacitive. Inter-dot 
charging energies (denoted by V) are assumed to be the same: they 
are characterized by the capacitance C, V = e /2C. Dashed lines 
indicate the electron transport through each dot. 

We model the TQD system (see Fig. 1) using a Hamiltonian 
consisting of three copies of the non-interacting resonant-level 



3 



model (each describing one QD and the effective single chan- 
nel of the electrons that the dot hybridizes with) and a cou- 
pling term which includes the inter-dot interactions and any 
possible inter-dot tunneling: 



H — Hi + H\ nt , 



(3) 



with 



and 



Hint = Wij n i n 2 + ^3 + nX 

(id) 



(5) 



Here cj. i is the creation operator for an electron with mo- 
mentum k in channel i, while d\ is the creation operator for 
an electron in dot i; the occupancy operator is defined as 
rii = d\di. Assuming flat conduction bands, the hybridisa- 
tion of each channel to the attached dot is characterized by a 
single number, I\ = npvf, where p is the density of states 
in the band which we take to be constant and of width 2D 
(flat-band approximation). Hereafter, we consider all energies 
in units of the half -bandwidth, D = 1. Vij is the charge re- 
pulsion between two dots, while tij is the hopping amplitude 
between two dots. For symmetrical configurations we sim- 
plify the notation as e = e^, t = t^, V = Vij. Notice that 
this model is similar to the Coqblin-Schrieffer SU(A) model. 
The Coqblin-Schrieffer model 54 describes an impurity in the 
N = 2j + 1 representation of the SU(7V) total angular mo- 
mentum group. The equivalence is established on low tem- 
perature scales where the charge fluctuations are quenched. 
Assuming that the conduction bands are particle-hole (p-h) 
symmetric (e_^ = — e^), the p-h transformation (d\ — > di, 



di -> d\, c\ 
constants) 



C-fe 



Ck — > c_ fc , etc.) leads to (up to irrelevant 



Hi =^2e k cl !i c kt i+(-e i )dldi + (-v i )Y^ (4,i^+ h -' 



and 



H^ = X) - - rij) - U 3 (d^ + h.i 



(id) 



Therefore 



(hi) 



(6) 



(7) 



(8) 



The model is p-h symmetric only for tij = 0, since finite 
inter-dot hopping breaks the bipartiteness. (The sign change 
of the hybridization Vi is of no physical consequence.) For a 
symmetric configuration, the model is p-h symmetric around 
the point e = — V. Therefore the parameter 6 = e + V is a 



measure of the departure from the p-h symmetry. It must be 
emphasized that the SU(3) Kondo effects that are discussed in 
the following occur away from the p-h symmetric point where 
the impurity charge fluctuates. This implies that the Kondo 
peak itself is not symmetric (confirmed by numerics, see be- 
low). 

The standard (Gell-Mann) parametrisation for the genera- 
tors of the SU(3) Lie algebra is K a = A a /2 with A a being 
Gell-Mann matrices: 



Hi = 53 e k4,i c M + e i4 d i + v i E \ c k,i di + hx -) ' ( 4 ) Ai 
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We thus define the SU(3) operators for the TQD system as 

ij 
ij 

O tota 1 = i m p + ^ ( fc)j 
k 

where a — 1 , . . . , 8, while i and j range over the three chan- 
nels, and k ranges over all conduction-band momenta. The 
Casimir operator of SU(3) is defined as 



KL,i = E(°a° ta1 ) 2 



(ID 



In a fully SU(3) symmetric case, the traces Tr(O^) are all 
equivalent. In numerical calculations, it is thus sufficient 
to calculate the expectation value of a single 0\ operator; 
the most convenient choice is 0\. The expectation value of 
(-^?otai) i s tnen 8 times this value. 

The behavior of an impurity system can be analyzed by 
studying its thermodynamic properties. In the following 
section we will consider the impurity orbital susceptibility 
Ximp(T 1 ) and the impurity entropy Si mp (T). These two quan- 
tities serve to establish the range of parameters for which the 
SU(3) spin Kondo physics is encountered. In the fundamen- 
tal representation of SU(3) one has (K 2 ) = 4/3. In the high- 
temperature regime where all eight dot states are equally prob- 
able, one has (K 2 ) = (6 x 4/3 + 2 x 0)/8 = 1, since there 
are six singly occupied states (by either one electron or by one 
hole) and two states corresponding to totally empty and totally 
full system. The impurity SU(3) orbital susceptibility (more 
precisely, this is the impurity contribution to the total system 
orbital susceptibility) is defined as 



Ximp 

(T) 



/3((At tal )(T) 



(A t 2 otal ) (T)) , 



(12) 
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where the bracket with subscript denotes the result for the 
system without the dots (i.e., the Hamiltonian H consists only 
of the conduction bands). Here j3 = 1/feeT with k B the 
Boltzmann constant. The value of k B Tx(T) therefore indi- 
cates the presence of a finite effective orbital local moment on 
the TQD and it can be used to classify the fixed point o 55 ' 56 . 

The impurity entropy is a measure of the number of the 
effective degrees of freedom of the TQD at a given parameter 
configuration. It is defined through 



(E -F) (E- F) 



T 



T 



(13) 



where E = (H) = Tv[H exp{-H/k B T)} and F = 
-k B T\nTi[exp(-H/k B T)}. 

We also compute the dot spectral functions A(ui,T) and 
compute the differential conductance through each dot using 
the Meir-Wingreen formula asS 



G(T) = Go 



duj 



7rI\4(w,T)dw, (14) 



where Go = e 2 /h and / = [l+expfw/fc^T)] -1 is the Fermi- 
Dirac distribution function; the chemical potential has been 
fixed at zero energy. 

The calculations have been performed using the numerical 
renormalization group metho d 55 ' 56 ' 58 ' 59 as implemented in the 
"NRG Ljubljana" code. We have used the discretization pa- 
rameter A = 8 with the z-averaging over N z = 8 values. We 
have verified that such a large value of A still produces reliable 
results by performing a convergence study as a function of A 
down to A = 2. In the NRG truncation, we have kept states 
with energy up to lOwyv where ljn is the characteristic energy 
scale at the iV-th NRG step, or at most 6000 states. For calcu- 
lating the spectral functions, we have used the complete Fock 
space metho d 60 - 61 . Very recently, a study of the fully symmet- 
ric SU(3) model has been performed with an implementation 
of the NRG which can explicitly use the SU(3) symmetry of 
the model to simplify the calculations^. Here we only use 
the U(l) total-charge-conservation symmetry, thus the calcu- 
lations are significantly more time-demanding. However, our 
approach makes it possible to study the effects of the sym- 
metry breaking terms, which is important for physical realiza- 
tions of this model. 



III. NUMERICAL RESULTS : VALENCE FLUCTUATING 
AND THE SU(3) KONDO REGIMES 

In order to identify the different regimes of the TQD system 
we consider the thermodynamic and transport properties. In 
Fig.|2]we show the basic results for the fully SU(3) symmetric 
case where all the dots and hybridizations are equivalent. The 
lead-dot couplings have common value of T = 0.01, and there 
is no inter-channel tunneling (i.e., t = 0). Since the system is 
symmetric with respect to the point e = — V for this choice of 
parameters, we consider only the value of the on-site energy 
e below —V (i.e., S < 0); other results can be obtained by an 
appropriate p-h transformation. 
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Figure 2: Temperature dependence of the impurity entropy (left pan- 
els) and the zero-temperature dot spectral function (right panels) for 
a range of the on-site energies e. We consider a symmetric triple 
quantum dot system. The total occupancy of the triple quantum dot, 
(n), and the zero-temperature linear conductance through one dot, 
G, are also shown. The interdot tunnelling is zero, t — 0. Other 
parameters are T = 0.01, V = 0.2. 



We observe that for a range of on-site energies —0.35 < 
e < —0.25, the occupancy reaches values close to 2. In 
this case, the system evolves from the free-orbital fixed point 
(fp) with impurity entropy of 3 In 2 to a local-moment fp with 
triple degeneracy (indicated by a In 3 plateau in the impu- 
rity entropy) as the temperature decreases below the charge- 
fluctuation scale of V. The triple degeneracy is finally 
lifted at low temperatures and then we reach a non-degenerate 
strong-coupling fp corresponding to the SU(3) Kondo regime. 
In this regime the zero-temperature linear conductance is near 
3/4, as predicted in Ref. [39| based on the Friedel sum rule 
arguments [see Eq. (1)]. 

A very different behavior is found near the p-h symmetric 
point at e = —V. Here the system evolves from the free - 
orbital to the valence-fluctuation fp with entropy In 6 (only 
visible as a weak bulge in the Si mp (T) curve in Fig. In 
this case the valence-fluctuation regime corresponds to charge 
fluctuations from n = 1 to n = 2 charge states. The entropy is 
eventually reduced from In 6 to zero at some low temperature. 
In this case the entropy is released as the system evolves to the 
strong-coupling fp without passing through the local-moment 
fp. 

By reducing the dot-lead hybridisation by half, i.e., for V = 
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Figure 3: (Color online) Temperature dependence of the impurity 
entropy (left panels) and the zero-temperature dot spectral function 
(right panels) for a range of the on-site energies e. Same parameters 
as in Fig. [2] but with smaller hybridization, V = 0.005. We plot the 
positive-frequency side of the spectral function on the logarithmic 
frequency scale. The dashed line corresponds to the characteristic 
conduction of G = (3/4) (e 2 /h). 



0.005, as shown in Fig. [3] the SU(3) Kondo regime is even 
more clearly discernible and we can see that the occupancy is 
pinned to the value 2 for a much broader range of dot potential 
energies. Here, the conductance reaches the universal value of 
G = 3/4(2e 2 /ft.) for a wide range of e due to a much more 
robust SU(3) Kondo state. 

In the SU(3) Kondo regime, the dot spectral density, shown 
in Fig.Hl displays a Kondo resonance with a maximum height 
shifted away from the Fermi level. As previously noticed 
this is the result of having the SU(3) Kondo regime away 
from the p-h symmetric point. In addition, the shifted spec- 
tral density produces a zero-temperature linear conductance 
G = 3/4(e 2 /h), which is by itself a hallmark of the occur- 
rence of the SU(3) Kondo physics. This is in contrast with the 
SU(4) case where the linear conductance coincides in value 
with the linear conductance for the SU(2) Kondo effect. 

The thermodynamic properties of the SU(3) Kondo regime 
are studied in more detail in Fig. [5] In order to gain some intu- 
ition about the role of the Kondo correlations in the SU(3) or- 
bital susceptibility we compare the case of an uncoupled TQD 
system, where all three lead-dot hybridizations are T = 0, 
with the case where the TQD is connected to leads and the 
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Figure 4: (Color online) Spectral function in the SU(3) Kondo 
regime. The triple quantum dot system is symmetric and the total 
occupancy is (n}=2. Note the pinning of the zero-frequency spectral 
density to the value A(0) = (3/4) /irT and the important fact that 
the actual Kondo peak is displaced away from a; = 0. 




Figure 5: (Color online) Thermodynamic properties in the SU(3) 
Kondo regime. We consider a symmetric TQD system. The total 
occupancy of the triple quantum dot, (n) = 2. We compare the 
cases of a decoupled triple quantum dot system (T — 0) and the 
triple quantum dot connected to the leads (T = 0.005) in which the 
Kondo correlations are present. 



6 



SU(3) Kondo state builds up. In the high-temperature limit, 
in both cases the spin susceptibility is 1 and all 8 TQD states 
are equally probable, thus there is In 8 = 3 In 2 impurity en- 
tropy. As the temperature decreases below V, so that the 
charge fluctuations are frozen out, the SU(3) susceptibility 
of 4/3 is established, as expected for the fundamental triplet 
representation of SU(3). This is the local-moment fp. The 
decoupled system stays in this fixed point down to T = 0, 
while at finite V the local moment is screened in the SU(3) 
Kondo effect and the susceptibility vanishes, as expected. 
At this point the system is in the non-degenerate strong- 
coupling fp in which the ground state corresponds to a SU(3) 
Kondo singlet state. Notice that the transition from the free- 
orbital to the local-moment regime and then eventually to the 
strong-coupling regime is fully analogous to the behaviour in 
the standard single-impurity Anderson model with the SU(2) 



symmetr 



,55,56.58 



We also emphasize that the low-temperature 



parts of the impurity susceptibility and impurity entropy are 
universal and that the scaling of the results for different pa- 
rameters is observed if the temperature axis is rescaled by an 
appropriately defined Kondo temperature Tk (see below). 

For completeness we also analyse the p-h symmetric point, 
in which the only low-temperature scale is T itself and there 
is no Kondo-like screening. In this model, the p-h symmet- 
ric point corresponds to a valence-fluctuation regime where 
charge fluctuations occur. In Fig. [6] the temperature depen- 
dence of the entropy is shown for a symmetric TQD and var- 
ious r values when e = — V. At high temperatures the TQD 
is found in the free-orbital regime where the TQD entropy is 
In 8. Then, the system crosses over on the temperature scale 
of V to a valence-fluctuation fp with a six-fold degenerate 
ground state in which the entropy reaches the value of In 6. 
In this case, there can be either a single electron or a single 
hole in the three dots for a total of six states with the same en- 
ergy. Decreasing further the temperature the system crosses 
over to the non-degenerate ground state with zero entropy at 
the temperature scale of T, see Fig. [6] There is no further 
dynamically generated low-energy scale in this case. 
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Figure 6: (Color online) Temperature dependence of the entropy at 
the particle-hole symmetric point for a range of hybridization param- 
eters F. We consider a symmetric triple quantum dot system with 
parameters V = 0.2, e = —0.3. 



The zero-temperature fixed points for different e form a 
line of fixed points which are related by the different strength 
of the residual potential scattering experienced by the quasi- 
particles. For an overview of the system behavior, in Fig.|7]we 
plot the zero-temperature total TQD occupancy, linear con- 
ductance, and charge fluctuations as a function of the on-site 
energy e for several choices of the hybridization T, The emer- 
gence of the Kondo plateau for low enough V is clearly visible; 
it coincides with the regions of low charge fluctuations in the 
TQD. 
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Figure 7: (Color online) Occupancy n, differential conductance G 
and charge fluctuations Sn 2 = (n 2 ) — (n) 2 as a function of the on- 
site energy for a range of hybridizations F. 



IV. EVIDENCE OF SU(3) KONDO CORRELATIONS IN 
TRANSPORT MEASUREMENTS 



In the previous section, we have demonstrated the occur- 
rence of the SU(3) Kondo effect by considering the thermo- 
dynamic properties. Usually, one way to probe the existence 
of Kondo correlations in QD systems is to measure the expo- 
nential dependence of the Kondo energy scale (ksTx) with 
the inverse hybridization 1/T as 



In TV oc 

r 



(15) 



In the NRG calculations this dependence is demonstrated by 
using a Tk defined from the entropy curve as S , i mp (Tx) = 



7 



O.lfcs, see the upper panel of Fig. [8] Plotting Tk for different 
r we uncover the exponential dependence, see the lower panel 
of Fig. m 
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Figure 8: (Color online) Top panel: Temperature dependence of the 
entropy for a symmetric triple quantum dot system for a range of 
hybridisation parameters Y. Other parameters are tuned so that the 
system is in the SU(3) Kondo regime, V = 0.2, e = —0.3. Bot- 
tom panel: Relation between the Kondo temperature Tk and the 
hybridisation T. Here we use an arbitrary definition of the Kondo 
temperature, Simp (Tic) = O.lfes. 

In view of this result, the occurrence of SU(3) Kondo cor- 
relations in a TQD experiment could be demonstrated by per- 
forming transport measurements, for instance by measuring 
the Kondo temperature Tk- Usually this is achieved by iden- 
tifying the half width at the half maximum (HWHM) of the 
non-linear conductance peak through one of the dots with Tk ■ 
This measurement would, however, be rendered problematic 
in this system due to the asymmetric shape of the Kondo res- 
onance. In dl/dV measurement with finite bias, the current 
is namely given approximately by 



I(T, V) = / [/lM - f R (w)]nrA(u, T)dw, 



(16) 



where we have neglected the voltage-dependence of the spec- 
tral function, and we assume = f(uj — V/2) and 
fu = f( u + V/2). In the zero-temperature limit, we thus 
find approximately 

dl r r A{V/2^)+A{-V/2,Q) 

av =Go7rT 2 ■ (17) 



Since the spectral function A(u>) is asymmetric, care is needed 
in extracting the width of the spectral function peak from 
the differential conductance measurements at finite bias volt- 
age (even when the non-equilibrium effects are neglected). 
When properly done, by measuring Tk and then varying T 
one should be able to get a relation that matches Eq. (TT3T >. 

Another signature of the SU(3) Kondo physics could be de- 
tected through the temperature dependence of the linear con- 
ductance G. Close to zero temperature, G is 3/4 in units of 
the conductance quantum Go = e 2 /h (note the absence of 
factor 2 in this spinless case). At small but finite temperatures, 
we find that the low-temperature conductance fits, similarly to 
the SU(2) Kondo case, to the empirical formula of the general 
form^ 



G(T) = G (l + (2 1 /* -\){T/T K f y 



(18) 



with s = 0.23 for the SU(2) case and with s = 0.28 for the 
SU(3) Kondo case. The value of s = 0.28 has been extracted 
from the NRG results for G(T). The fitting to G(T) is per- 
formed for a symmetric TQD at e = —0.2, and T = 0.005, 
see Fig. [9] It is interesting to notice that a single fit for- 
mula, Eq. ( fT8l . applies over many orders of magnitude in 
temperature in a number of quantum impurity models that ex- 
hibit Kondo effects of very different kinds; see, for example, 
Refs.iHi 
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Figure 9: (Color online) Temperature dependence of the linear con- 
ductance in the SU(3) Kondo regime of a symmetric triple quantum 
dot. 



DEPARTURE FROM THE SU(3) STRONG-COUPLING 
FIXED POINT 



In real experiments it is extremely challenging to construct 
perfectly symmetric multi-dot systems. Therefore, in order 
to experimentally detect SU(3) Kondo correlations we need 
to know to what extent the SU(3) Kondo physics is robust 
against all possible local perturbations. These perturbation 
can arise from asymmetric lead-dot tunneling couplings, pos- 
sible tunneling events among the dots, and different on-site 
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dot level potentials or distinct inter-dot Coulomb energies. In 
spite of the presence of unavoidable local perturbations in real 
setups that eventually destroy the SU(3) Kondo state, below 
we give a protocol to restore the Kondo correlations by prop- 
erly adjusting the dot potentials. 



A. Asymmetric lead-dot hybridizations 

First, we analyse the effect of having asymmetric lead-dot 
coupling. If one of the hybridizations Ti, say T\, is made 
weaker, the SU(3) symmetry is broken. Even small changes 
of Ti are sufficient; the effect is similar to the induced mag- 
netization by a ferromagnetic conduction band in the spinfull 
single QD device. In the upper panel of Fig.[lO]we show the 
impurity entropy evolution as the temperature is lowered for 
constant r 2 = r 3 = T and Ti ranging from I 1 ! = T (sym- 
metric configuration) to Ti = 0.9T (asymmetric lead-dot cou- 
plings). When the temperature is lowered, in the asymmetric 
lead-dot configuration the system flows from the SU(3) local- 
moment fp (with impurity entropy In 3) to a new SU(2) local- 
moment fp with two-fold degeneracy (with impurity entropy 
In 2). The SU(2) local moment is then screened in the con- 
ventional SU(2) Kondo effect which lifts the degeneracy at 
temperatures well below the new SU(2) Kondo energy scale 
(see lower panel in Fig. [TU]for the comparison between the 
entropy evolution for the asymmetric coupled triple dot case 
and the universal SU(2) Kondo model). If, however, Ti is 
increased above T 2 = T 3 rather than decreased, the system 
crosses over from the SU(3) local-moment fp to the frozen- 
impurity fp without any Kondo screening (results not shown). 
The same happens when all three Ti are different. 

As already mentioned, the asymmetry of the lead-dot cou- 
plings in the SU(3) Kondo effect is analogous to the SU(2) 
Kondo physics in the presence of ferromagnetic contacts. In 
the latter case, the Kondo resonance is split due to the ap- 
pearance of an induced exchange field because of the polar- 
ized contacts. The same physical behavior is obtained for the 
SU(3) Kondo case. Figure QT] shows this result. We plot the 
spectral densities for the dot 1 and dot 2, denoted as A\ (ui) 
and A 2 (w). We consider the case of symmetric couplings 
(r = Ti = T 2 = r 3 ) in which Ai(ui) = A 2 (u) and the 
case of asymmetrical lead-dot coupling configuration where 
r 2 = T 3 = 0.005 and Ti = 0.003. In the latter case the two 
spectral functions show SU(3) Kondo-peak splitting which is 
better visible in the close-up shown as an inset in Fig.QT] It is 
noteworthy that the Kondo spectral peaks for A 2 = A3 reach 
a high value (approaching, in fact, the unitary limit), while 
that for Ai is strongly suppressed. This is related to the fact 
that the dots 2 and 3 are SU(2) Kondo screened, thus their 
zero-bias conductance remains nearly unitary, while the dot 1 
becomes decoupled and it is only weakly conducting. 




Figure 10: (Color online) Top panel: Temperature dependence of the 
impurity entropy for a range of the hybridization parameter Ti, while 
r2 = T3 — r is held fixed. Bottom panel: Impurity entropy for a 
strong symmetry breaking and a fit of the SU(2) Kondo screening 
cross-over with the universal SU(2) entropy curve for the spinfull 
single dot case. 
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B. Unequal on-site dot energies 



The SU(3) Kondo state can also be destroyed by having 
non-equal dot level energies. This case is illustrated in Fig. [12] 



Figure 1 1 : (Color online) Spectral functions in the case of non equal 
hybridisations Ti . The dashed line corresponds to the fully symmet- 
ric case with T = T t = 0.005. 
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where it is shown how the SU(3) symmetry is broken by 
changing the on-site energy e\ away from the common value 
e, i.e., ei = e + 5ei, with 5ei being the detuning. The result- 
ing new state depends on the direction of the detuning. For 
positive detuning 5e\ > 0, shown in the upper panel in Fig. 
[T2l the SU(3) Kondo effect is quenched on the energy scale 
of the detuning. However, for negative detuning 8e\ < 0, 
shown in the lower panel, the system evolves from the three 
fold degenerate local moment fp with 5i mp = In 3 to a two 
fold degenerate local moment fp in which 5j mp = In 2. In 
this case the SU(2) local moment emerges from two states 
that originally formed the SU(3) triplet local moment. 




Figure 12: (Color online) Temperature dependence of the impurity 
entropy for a range of the parameter 8e\ defined as ei = e + <5ei, 
while £2 = £3 = e with e = —0.3. The rest of the parameters corre- 
spond to a symmetric triple quantum dot configuration with V = 0.2, 
T = 0.005. 



C. Finite interdot tunneling: t / 

We have also investigated the fact that the three-fold sym- 
metry can also be broken by a finite interdot tunneling. The 
behavior yet again depends on the sign of t: the ground state 
may either be a single state (not Kondo screened) or a two-fold 
degenerate pair which is Kondo screened. The positive t case 
is considered in the upper panel of Fig. Qj] Here the entropy 
evolves from In 3 to zero on the temperature scale of t signal- 
ing the destruction of the SU(3) Kondo singlet. For negative 
t, the lower panel in the figure shows that the impurity entropy 



evolves from In 3 to In 2 on the temperature scale of \t\ and the 
latter corresponds to local moment fp of two-fold pair of states 
and it constitutes a SU(2) local-moment which undergoes the 
SU(2) Kondo screening at much lower temperatures. 
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Figure 13: (Color online) Temperature dependence of the triple 
quantum dot entropy for different values of inter-dot tunneling. 



D. Restoration of the SU(3) Kondo physics 

Symmetry can also be broken by having different inter-dot 
charge repulsion parameters Vij . In experiments, the capac- 
itive couplings Vij are the most difficult to control, followed 
by the inter-dot hopping parameters and hybridizations I\, 
while the on-site energies e; are typically the easiest to tune. 
There are eight "directions" for an SU(3) symmetry-breaking 
field, corresponding to the eight generators of this symme- 
try. A consideration of the Gell-Mann matrices in Eq. (O 
suggest that two (A3 and As) are associated with the energy 
levels, since they are diagonal, while the remaining six are 
out-of-diagonal and thus associated with the inter-dot tunnel- 
ing. This immediately suggests that symmetry breaking by 
tunneling cannot be compensated electrostatically. It is thus 
essential, first of all, to ensure that the inter-dot tunneling is 
low-enough for the observation of the SU(3) Kondo effect. 
In essence, the tunneling rate should be much smaller than 
the anticipated Kondo temperature scale (in appropriate fre- 
quency units). Any remaining asymmetry then corresponds 
to SU(3) "fields" in A3 and Ag directions which arise from 
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asymmetric ej, I\ and Vij. Assuming that the three on-site 
energies can be freely and independently tuned, it appears 
possible to compensate the asymmetries in I\ and Vij since 
there are three parameters to drive two "fields" to zero. We 
demonstrate this procedure in Fig. [14] where an asymmetry 
in the hybridization constants is compensated by tuning the 
on-site potentials. Thus, making use of the high degree of 
tunability in QD devices, any source of symmetry breaking 
which naturally arise from the manufacturing process can be 
compensated by properly adjusting the gate voltages. 




Figure 14: (Color online) Restoration of the SU(3) Kondo effect in 
an asymmetrically lead-dot coupled triple quantum dot (Ti = 0.004, 
and r2 = T3 = 0.005) by tuning the on-site dot potential ei = 
e + 5ei, with V = 0.2, and e — —0.3. The SU(3) symmetry is 
effectively restored for 5ei = 3.25 10~ 5 , which is indicated by the 
impurity curve having the universal SU(3) shape. The curves for 
parameters near the restoration point show some wiggles; these are 
numerical artifacts. 



VI. SPINFULL TRIPLE QUANTUM DOT 

For completeness, in this section we briefly discuss the ver- 
sion of our model with spin degrees of freedom and additional 
on-site electron-electron repulsion terms U . The Hamiltonian 
takes the form of three copies of the single-impurity Anderson 
model and a coupling term with all inter-dot terms: 



(19) 



i=l 



with 



k.rj 



V i E ( 4,<r,i d <M + nC 



(20) 



and 



flint = E 



k.a 



1 '<./''<»./ + Uj E \ d l,i d °,i + h c 



(21) 



Most terms have the same meaning as in the spinless model, 
Eq. (O, but now n a ,i = d a and rii = n-j-^ + n^i. 

Here we will only consider some features of this model. We 
restrict our attention to a model with no inter-dot tunneling, 
tij = 0, and full orbital symmetry, i.e., = e, Ui = U, 
= V and = v (or, equivalently, = T). Despite 
the high symmetry of this problem, it is still too complex for 
a detailed study using the NRG technique. For this reason, 
we resort to a different numerical approach, the Gunnarson- 
Schonhammer variational method. Following Refs.lr34landl65l 
we form a variational ansatz for the ground state wavefunction 
|0) of the Hamiltonian ( fl"9l l. 



|0) 



+ 



+ 



^ . ^riin2Ti3-fniri2n3 |0) ~l~ 

3 

a 

4. 4 c ^|6)-(22) 



E 
E 



/V nin 2 n 3 - r ™i™27i3 / t ' 



Here 1 0) is the ground state wavefunction of the noninteract- 
ing part (U = V = 0) of the Hamiltonian (fl~9b with renormal- 
ized dot energy levels e and lead-dot couplings v, 



H(e,v) 



E 



+ E 

i 

h.c 



EH 



j c ka,' 



(23) 



Projectors P ni n 2 n 3 project this state to subspaces with rii = 
0,1,2 electrons in the i-th dot. Variational terms in the sec- 
ond and the third row of Eq. (1221) provide states containing an 
electron above the Fermi energy and a hole below the Fermi 
energy in one of the leads, respectively. An approximation to 
the true Hamiltonian ( fl~9l ground state energy Eq (e, v) and 
the corresponding ground state wavefunction coefficients A 
are calculated by solving the Schrodinger equation within the 
Hilbert space of the ansatz. The ground state energy is then 
further minimized with respect to e — > eo an d v — > ^0* provid- 
ing us with the noninteracting part of the Fermi liquid quasi- 
particle Hamiltonian H (eo, vq) of our problem which we use 
to calculate the zero temperature conductance. The occupancy 
of the dots and its fluctuations are calculated from the ground 
state wavefunction |0). 

To study the interplay of the inter-site and on-site charge 
repulsion terms U and V, we first fix U and increase V. For 
V = 0, the system consists simply of three copies of the 
single-impurity Anderson model, the properties of which are 
well known^S. In the interval —U + T < e < -T, we 
expect the emergence of the SU(2) Kondo effect in each of 
the three channels independently, thus a plateau of unitary 
conductance through each dot (note that the unitary limit is 
now 2e 2 /h due to the spin factor). For non-zero but moder- 
ate V < U, we expect the occurrence of collective Kondo 
screening which affects all three quantum dots. In this case, 
in addition to the spin degree of freedom on each dot, there 
is an orbital degree of freedom, as in the spinless case that 
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has been discussed in the previous sections. The two de- 
grees of freedom may become intertwined like in the SU(4) 
Kondo effect in carbon nanotubes where the spin and isospin 
degrees of freedom are combined. Here, however, the sym- 
metry is SU(2) sp in x SU(3) or bitai and richer behavior is ex- 
pected. This model is closely related to the impurity mod- 
els studied in the context of dynamical mean-field theory for 
correlated bulk metals, in particular for transition-metal com- 
pounds where the three-fold degenerate t<z g d-electron orbitals 
play the main role. The results in Fig. [T31 indicate the emer- 
gence of three different kinds of the Kondo plateau for on- 
site energies corresponding to the total TQD occupancy be- 
ing pinned to an integer value. At the p-h symmetric point, 
we observe a Kondo plateau for all values of V < U; this 
is thus contrary to the behavior found in the spinless model 
where in this regimes the valence fluctuates and the Kondo 
effect does not occur. Two additional types of the Kondo 
effect are present, one for single-electron (single-hole) occu- 
pancy, (n) = 1, 5, characterized by the conductance pinned to 
G = 2e 2 //isin 2 [(7r/2)(l/3)] = 0.5(e 2 /h), another for two- 
electron (two-hole) occupancy, (n) =2,4, characterized by 
G = 2e 2 //isin 2 [(7r/2)(2/3)] = 1.5(e 2 //i). In the first case, 
the single electron has both spin and orbital degree of freedom 
(both in the fundamental SU(2) and SU(3) representations), 
thus this corresponds to the SU(6) Kondo effect. In the second 
case, two electrons are comined into a more complex object 
which is then Kondo screened. 



VII. CONCLUSIONS 



We have demonstrated the formation of the SU(3) singlet 
Kondo state in a spinless triple quantum dot system. We 
have performed a thermodynamics analysis with the help of 
the magnetic susceptibility (ximp) and the impurity entropy 
(Simp)- For a highly symmetric triple quantum dot we find, 
when the dot system charge is either n = 2 or n = 1 (away 
from the particle-hole symmetric point for which n = 1.5), 
that the system evolves from the free-orbital regime at high 
temperatures with Tximp = 1 and Si mp = In 8 to the three- 
fold degenerate local moment fixed point with Tximp = 4/3 
and Si m p — In 3 towards a non-degenerate strong coupling 
fixed point where the SU(3) spin is totally screened. This 
cross-over occurs on an exponentially low temperature scale, 
IuTk oc — 1/r. In contrast, for the electron-hole sym- 
metric point, the system evolves from the free-local moment 
regime to a valence-fluctuating fixed point where the entropy 
Si m p = In 6, and then to zero on the scale of T. Addition- 
ally, we have investigated possible perturbations that affect 
the formation of the SU(3) singlet Kondo state in real se- 
tups. Among these perturbations we have studied how non- 
symmetrical lead-dot couplings or on-site potential energies, 
and finite inter-dot tunneling rates break the SU(3) Kondo 
physics. We have described the procedure to restore the SU(3) 
Kondo physics in real systems by properly tuning the on-site 
dot energies. Finally, we described the more complex behav- 
ior of the spinfull case where the Kondo effect occurs for all 
integer occupancies of the triple quantum dot, including at the 
particle-hole symmetric point. 




Figure 15: (Color online) Occupancy n, differential conductance G 
and charge fluctuations Sn = {n 2 } — {n) 2 as a function of the on- 
site energy e = — U /2 — 2V + 5 for on-site repulsion U — 0.32 and 
for a range of interdot repulsions V. V = 0.01. 
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